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ABSTRACT 

The presence of an extended blue horizontal branch (HB) in a stellar population is 
known to affect the age inferred from spectral fitting to stellar population synthesis 
models. This is due to the hot blue component which increases the strength of the 
Balmer lines and can make an old population look spuriously young. However, most 
population synthesis models still rely on theoretical isochrones which do not include 
realistic modelling of extended HBs. In this work, we create detailed models for a range 
of old simple stellar populations (SSPs), with metallicities ranging from [Fe/H]=— 1.3 
to solar, to create a variety of realistic HB morphologies, from extended red clumps, to 
extreme blue HBs. We achieve this by utilising stellar tracks from the BaSTI database 
and implementing a different mass loss prescription for each SSP created. This in- 
cludes setting an average mass and a Gaussian spread in masses of individual stars 
coming on to the zero age HB for each model, and hence resulting in different HB 
morphologies. We find that, for each metallicity, there is some HB morphology which 
maximises H/3, making an underlying 14 Gyr population look ^5 — 6 Gyr old for the 
low and intermediate metallicity cases, and as young as 2 Gyr in the case of the solar 
metallicity SSP. We explore whether there are any spectral indices capable of break- 
ing the degeneracy between an old SSP with extended blue HB an d a truly young or 
intermediate age SSP, and find that the Can index of iRosd (|1984f) and the strength 
of the Mgii doublet at 2800A are promising candidates, in combination with H/3 and 
other metallicity indicators such as Mg6 and Fe5406. We also run Monte Carlo simula- 
tions to investigate the level of statistical fluctuations in the spectra of typical stellar 
clusters. We find that fluctuations in spectral indices are significant even for average 
to large globular clusters, and that various spectral indices are affected in different 
ways, which has implications for full-spectrum fitting methods. Hence we urge cau- 
tion if these types of stellar clusters are to be used as empirical calibrating objects for 
various aspects of SPS models. 

Key words: stars: horizontal branch - globular clusters: general - galaxies: stellar 
content. 



1 INTRODUCTION 

In recent years, stellar population synthesis (SPS) mod- 
els have rapidly become a fundamental tool in the 
study of both Galactic an d extr a galactic stella r pop- 
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ulations (see Percival et all 20091; 

Coelho et al. I 120071: Vazquez fc Le ithcrcr 2005; Mar astonl 
20051 : iBruzual fc Charlotl 120031 for some recent examples) . 
Despite the increasing sophistication of the underlying stel- 
lar models used in SPS and ever-growing libraries of spec- 
tra (both empirical and synthetic) available to modellers, 
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there are still areas in which models, and methods for fitting 
them to observational data, can be significantly improved 
and expanded. Some of the current shortcomings are due to 
the need to adopt simplifying assumptions about the stel- 
lar population in question in order to make some problems 
tractable, e.g. the fitting of simple stellar population (SSP) 
models to galaxies, assuming that they can be represented 
by a single age, single metallicity model. 

A key area in which simplifying assumptions are im- 
plicitly made is the morphology of the horizontal branch 
(HB) for individual SSPs - this is largely due to the limita- 
tions of theoretical isochrones, which are used as the basis of 
all SPS models (except those of Maraston and co-workers, 
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who use a fuel-consumption based method - see iMarastonl 
120051 and references therein) . However it is well known that 
the morphology of the HB can impact strongly on the inte- 
grated light of stellar populations especially if an extended 
blue component is present, potentially affecting both colours 
and line indices, and he nce impacting on inferred ages for 
these systems (see e.g. Ocvirkl 2010l; Conrov et al.l 120091 ; 
ISchiavon et~ai1l2004l ; iLee et al.ll2002l . l2000h . Hence it is vital 



to assess the ways which these 'simplified' HBs inherent in 
theoretical isochrones impact on SPS models, and whether 
a more realistic and detailed treatment of the HB is nec- 
essary, or indeed practical. This is the purpose of the work 
presented here. 

For unresolved stellar populations our knowledge of 
their ages and elemental abundances is largely derived 
through the fitting of diagnostic spectral indices, hence the 
focus of our work is to investigate the effect of a detailed 
treatment of the HB on key spectral indices, such as H/3. 
However, it is important to realise that our results will also 
impact on analyses which use full-spectrum fitting methods, 
as will be discussed in Section 

In real stellar systems the horizontal extension of the 
HB is governed by stochastic mass loss in stars approaching 
the tip of the red giant branch (TRGB). All the stars in a 
particular SSP leave the TRGB with the same core mass 
but individual stars have different masses remaining in the 
outer layers, which determine the position of the star along 
the zero-age horizontal branch (ZAHB). In stellar models 
this mass loss is parameterised by the mass loss parameter, 
77, which is assigned some fixed v alue calculated according 
to the Reimers mass loss relation (|Reimerslll975h . Hence in 
theoretical isochrones the HB comprises a single mass point, 
with no extension. 

It is worth remembering here that an isochrone con- 
sists of a series of evolutionary points (EPs) which define 
the locus of points for an SSP. Whilst each EP defines the 
appropriate stellar parameters for any star located at that 
point, in terms of mass, effective temperature and luminos- 
ity, they do not represent individual stars themselves. In or- 
der to create an SSP, an isochrone is "populated" according 
to some initial mass function (IMF) which effectively gives 
the appropriate weighting to each EP. In most population 
synthesis work the isochrone is treated as an analytical func- 
tion for this purpose, so that all EPs along the isochrone are 
smoothly populated and the ZAHB remains as a single mass 
point. In order to create a model with an extended HB, effec- 
tively incorporating a spread in mass loss, the isochrone (or, 
at least, the HB portion) must be populated with a discrete 
number of stars, so that each star in the HB phase can be 
assigned a specific mass from within some range of masses 
- this requires interpolation between individual core-helium 
burning stellar tracks. 

Creating an integrated spectrum for these extended 
HB models is a computationally expensive procedure since 
the time taken is proportional to the number of points 
in the simulation (for details on how integrated spectra 
for SSPs are produ ced the interested reader is referred to 
iPercival et al. 1 120091 . hereafter P09). For the analytical case, 
the number of points is just the number of EPs along the 
isochrone, which varies between a few hundred and around 
two thousand depending on which isochrone set is used. To 
populate an isochrone with individual stars, one has to con- 



sider the problem of statistical fluctuations which are likely 
to arise as a result of low numbers of stars (and hence poor 
sampling) in the later stages of evolution, including the HB 
phase. This is of particular relevance in the case of an SSP 
with an extended blue HB, since these hot blue HB stars give 
rise to strong Balmer lines, which are generally used as the 
primary age indicators for stellar populations. This problem 
can only be overcome by using large numbers of points in 
the simulation, potentially up to ~ 10 . In fact, part of the 
work presented here is to explore how large these statistical 
fluctuations can be in real stellar clusters, and how many 
input stars are needed in the models to avoid significant un- 
certainties in ages and metallicities derived from the final 
integrated SSP spectra. 

Another problem in creating models with extended HBs 
is that there is still no theory which predicts mass loss rates 
from stellar parameters and so the value of 77 is chosen ar- 
bitrarily, usually to reflect the typical HBs seen in Galactic 
globular clusters. This is complicated by the fact that clus- 
ters apparently with the same age and chemical composition 
can have different HB morphologies. However, in principle 
77 can take any value and so producing a database of SPS 
models to cover all possibilities is not fea s ible. 

The pioneering work of ILee et all <|2000f > (hereafter 
L2000) highlighted the importance of including realistic HBs 
in stellar population work. L2000 created models for 15 SSPs 
using a single value for the mass loss parameter for all their 
models, which was chosen to replicate the range of HBs 
seen in Galactic globular clusters. They included a Gaus- 
sian spread in the HB mass distribution of <jm=0.02Mq to 
simulate the observed extension of the HB. This work graph- 
ically demonstrated the systematic variation of HB morphol- 
ogy with age and metallicity (see their Figure 5), showing 
that, at fixed 77, HBs generally become bluer as metallicity 
decreases and as age increases. L2000 note that the strength 
of the H/3 line in the integrated population does not in- 
crease monotonically as metallicity decreases at a fixed age, 
but peaks at some intermediate metallicity and then falls 
again as metallicity continues to decrease. This is because 
H/3 reaches a maximum strength in stellar spectra when the 
effective temperature, T e //, is around 9500K and so H/3 is 
maximised if the distribution of stars along the HB centres 
around this value. In these cases the contribution to H/3 
from the HB completely dominates over the contribution 
from the main-sequence turnoff and makes the population 
look spuriously young. 

It is apparent from the L2000 figures, but not specif- 
ically noted by them, that the actual extension (i.e. from 
red to blue) of the HB in the colour-magnitude diagrams is 
very different for the various SSPs, even though the value 
of 77 and om are the same. This is because the stellar ef- 
fective temperature in this phase is extremely sensitive to 
very small changes in envelope mass and so the resulting 
HB morphologies can be very different, even when the same 
mass loss prescription is used. 

In this paper we use similar techniques to L2000 to 
model extended HBs, populating them with discrete num- 
bers of stars by interpolating between core-helium burning 
stellar tracks and using a Gaussian spread in the mass loss 
distribution, and we extend that work in 3 key ways. Firstly 
we explore the effects of varying rj at fixed age and metal- 
licity, which enables us to model high mass loss in metal 
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Figure 1. 4 models from the ct-enhanced, [Fe/H]=— 0.70 set. For 
clarity, only 500 points are plotted for each HB (~10000 were used 
in each simulation). The 14 Gyr isochrone used for the underly- 
ing population up to the TRGB is also plotted. For reference 
purposes, the later evolutionary stages of the isochone are also 
shown for the 2 'standard' cases in the BaSTI database, i.e. with 
fixed r)=0.2 and r;=0.4. Horizontal branch morphology is param- 
eterised by HB type=(B - R)/(B + V + R) - see text for details. 



rich systems to see whether we can mimic a system such as 
M32, which has approximately solar metallicity, but an ex- 
tended blue HB component. We are also able to create SSPs 
at different metallicities with bimodal HBs, similar to those 
observed in several Galactic globular clusters. Secondly, we 
create several models which we populate with numbers of 
stars representative of typical stellar clusters, and use Monte 
Carlo techniques to explore the impact of statistical fluc- 
tuations on measured diagnostic line indices, such as H/3. 
Thirdly, we create full integrated high resolution spectra for 
all our models and simulations which enables us to assess 
the behaviour of all diagnostic line indices (within a wave- 
length range 2500 A to ~6000A), as well as the continuum 
flux. We also investigate whether there are any diagnostic 
indices capable of breaking the degeneracy between an old 
SSP with extended blue HB (hence strong H/3) and a truly 
young or intermediate age SSP. 



2 MODELLING SSPS WITH EXTENDED 
HORIZONTAL BRANCHES 

2.1 Method 

As a starting point, we created 16 SSPs with extended HBs, 
in 4 sets of 4 models each: these comprise 2 sets with scaled- 
solar abundance ratios, at [Fc/H]= —1.27, +0.06 (solar 
metallicity) and 2 sets with a-enhanced ratios, at [Fe/H] = 
— 1.31, —0.70. In each case the underlying population was 
created from an appropriate metallicity 14 Gyr isochrone 
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Table 1. Parameters used for a-enhanced HB models. 



Model 


[Fc/H] 


< M > (M ) 


CM 


1 


-1.31 


0.70 


0.02 


2 


-1.31 


0.64 


0.03 


3 


-1.31 


0.58 


0.03 


4 


-1.31 


0.51 


0.005 


1 


-0.70 


0.71 


0.02 


2 


-0.70 


0.60 


0.03 


3 


-0.70 


0.55 


0.02 


4 


-0.70 


0.495 


0.005 


ldD16 Zi. 


Parameters used for scaled-solar H 


Model 


[Fc/H] 


< M > (Mq) 




1 


-1.27 


0.72 


0.02 


2 


-1.27 


0.65 


0.02 


3 


-1.27 


0.58 


0.02 


1 


-1.27 


0.53 


0.01 


1 


+0.06 


0.70 


0.03 


2 


+0.06 


0.63 


0.03 


3 


+0.06 


0.56 


0.02 


4 


+0.06 


0.51 


0.01 



from the BaSTI databas^H ( Pietrinferni et al ] |2004l . I2006T I. 
For each model the isochrone was populated analytically up 
to the TRGB and the integrated spectrum for this portion 
of th e isochrone was cre ated using high resolution spectra 
from lMunari et al.l (120051 ) . as described in P09. The HB was 
treated separately and populated with a large number of in- 
dividual points (> 10000, in order to avoid the problem of 
statistical fluctuations) with some mean mass (correspond- 
ing to some value of rf), and with a Gaussian spread in the 
mass distribution, <tm, by interpolating between core-helium 
burning stellar tracks. Spectra were assigned to each point 
in the HB simulation as usual, i.e. matching T e //, log g, and 
[Fe/H], and these were then summed together. The summed 
HB spectrum was then scaled appropriately before being 
added to the spectrum of the underlying population to cre- 
ate the final integrated spectrum for each model. 

Equivalent widths of various diagnostic line indices were 
measured on the resulting integrated spectra, as described in 
P09, so that comparisons could be made with the BaSTI SPS 
database. Briefly, the line indices dis cussed here a r e from 
the Lick/IDS bandpasses tabulated in iTrager et all l|l998l) , 
unless otherwise stated. The quoted line strengths were ob- 
tained using the LECTOR programme by A. Vazdekifl and 
are those directly measured on the spectra, i.e. they are not 
transformed onto the Lick system. We focus in particular 
on the H/3 index in the following section since it is used 
primarily as an age indicator for SSPs. 

In each set of models, the 4 cases were created to ap- 
proximately match the extremes of HBs seen in Galactic 
globular clusters, from a purely red, but extended HB (i.e. 
a red clump - HB model 1) to an extremely blue HB with 
extended blue tail (HB model 4), plus 2 intermediate cases 

1 Available at |http: / / afbione.oa-teramo.inaf ,it/| 

2 See http://www.iac.es/galeria/vazdekis/index.html 
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Figure 2. H/3 vs. Fe5406 a-enhanced grid with results for the 
2 sets of a-enhanced models overplotted. The underlying grid 
is from the BaSTI SPS database, with fixed r) = 0.2, for SSP 
ages 1.25, 3, 6, 8, 10 and 14 Gyr (from top to bottom) and 
[Fe/H]=-1.84, -1.31, -1.01, -0.70, -0.29, +0.05 (increasing 
from left to right). The dashed line is for the 14 Gyr BaSTI SSPs 
with fixed rj = 0.4. Cross es and asteri s ks are observational data 
for Galactic GCs from ISchiavon et al. - the asterisks de- 

note clusters with HB type > 0.8 (i.e. predominantly blue). 



(HB models 2 and 3). We also combined spectra from the 
extreme red and extreme blue cases to create models with 
a bi- modal HB. Parameters used for each set of models, in 
terms of metallicity, mean mass < M >, and om, are tab- 
ulated in Tables 1 and 2, respectively, for the a-enhanced 
and scaled-solar sets. 

Figure [1] shows the V/ (B-V) colour magnitude diagrams 
(CMDs) for the a-enhanced, [Fe/H]= -0.70 set of 4 models, 
which illustrates the typical HB morphologies for each set. 
T he HB morpholog ical type (B-R)/(B+V+R), as defined 
m Lee et al l (| 19941 ) . was calcul ated using the boun daries of 
the instability strip defined in ICordier et al.l |2007l ). and is 
notated on each plot. The appropriate isochrone (used to 
create the underlying population up to the TRGB) is also 
plotted for the 2 'standard' cases in the BaSTI database, i.e. 
with fixed r\ of 0.2 and 0.4. At this intermediate metallicity 
the ZAHBs for the 77=0.2 and 77=0.4 cases have similar ef- 
fective temperature, hence there is only a small difference in 
their location in the CMD. The 77=0.4 case always produces 
a hotter (i.e. bluer) ZAHB then the t;=0.2 one, but for metal- 
licities of [Fe/H]~ —0.70 and greater, this still produces a 
ZAHB which is cooler (i.e. redder) than the turn off (TO), 
even at old ages. However, for metallicities lower than this, 
at old ages, the 77=0. 4 case produces a ZAHB which is very 
different from the 77=0. 2 case, and is much hotter (bluer) 
than the TO. The difference between the 2 cases, and the 
impact on derived ages, is discussed in Section \2. 2 1 




[MgFe] 



Figure 3. H/3 vs. [MgFe] scaled-solar grid with results for the 2 
sets of scaled-solar models overplotted. The underlying grid is as 
in Figure [2] but for [Fe/H]=-1.79, -1.49, -1.27, -0.96, -0.66, 
-0.35, +0.06, +0.40. 



2.2 Results 

Figure [2] shows results from the a-enhanced HB models in 
the H/3 v. Fe5406 diagnostic grid, i.e. where H/3 is the age in- 
dicator and Fe5406 traces iron abundance, [Fe/H]. Similarly, 
Figure [3] shows results from the scaled-solar HB models in 
the H/3 vs. [MgFe] grid, where [MgFe] is a tracer of total 
metallicity, Z (see P09). In these figures, HBl refers to the 
red clump only HB model, HB4 refers to the extended blue 
tail case and HB2 and HB3 are the intermediate cases. The 
underlying models from which these grids are constructed 
are the BaSTI 77=0.2 models (solid lines), with the 77=0. 4 
models at 14 Gyr only overlaid (dashed lines), described in 
P09. 

Before discussing the results of the extended HB simula- 
tions it is important to note the behaviour of H/3 in the fixed 
?7=0.4 models, at old ages, compared to the 77=0.2 ones. Fig- 
ures [2] and [3] clearly demonstrate that at metallicities above 
[Fe/H]~ —0.7, the H/3 line has similar strength for the 2 
cases, indicating that their HBs have similar average tem- 
peratures and both appear as red clumps in a CMD. How- 
ever between [Fe/H]= -0.7 and [Fe/H]= -1.0, the H/3 line 
suddenly becomes much stronger for the ?7=0.4 case, which 
is an indication that the HB stars are much hotter than for 
the 77=0.2 case, and hence the HB would appear much bluer 
than the turn off in a CMD. For the 77=0.4 models, H/3 peaks 
at a metallicity of around [Fe/H]= —1.0 and then decreases 
again as the metallicity continues to decrease. This is the 
same behaviour noted by L2000 and happens because the 
average temperature of the HB stars gets increasingly hot 
as the metallicity decreases until it is significantly hotter 
than 9500K (where the H/3 strength peaks). 

In all 4 model sets, an extended red only HB (model 
HBl) only negligibly affects H/3 compared to the 'standard' 
BaSTI SSP model with single mass ZAHB and 77=0.2, hence 
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the inferred age from these diagnostic grids is the correct 
one for the underlying stellar population, i.e. 14 Gyrs. For 
the 3 sets of models with sub-solar metallicity (i.e. the 2 a- 
enhanced sets and one of the scaled-solar sets), H/3 reaches 
a peak for one of the intermediate morphology cases (model 
HB3), and is similar to the fixed 77=0.4 case, whilst for the 
scaled-solar solar metallicity set, H/3 peaks for model HB4, 
i.e. the extended blue tail case. For the lower metallicity 
models, at [Fe/H]~ —1.3, the peak in H/3, whether from the 
extended HB model or from the fixed 77=0.4 case, implies 
an age around 5-6 Gyrs, whilst for the solar metallicity case 
the inferred age is <3 Gyr (remembering that the underlying 
population is a 14 Gyr one in all cases.) 

Looking in more detail at the individual models, it is 
evident that if a particular HB simulation has a distribution 
with a mean mass which corresponds to T e // ~9000-9500K, 
then H/3 is maximised, implying spuriously young ages from 
the integrated SSP spectra (confirming the result noted by 
iLee et allEiooh . 

We created bi-modal HBs by combining models HB1 
and HB4 in each set, and varying the fractions of the red and 
blue components. In all cases it was found that the change in 
H/3 strength just scales linearly with the change in blue HB 
fraction. This is illustrated in Figures[2]and[3]where the open 
star symbol represents the 50/50 case (i.e. 50% red, 50% blue 
HB) which always lies exactly half way between model HB1 
and model HB4. Therefore inferred ages for an underlying 
14 Gyr population can lie anywhere between <3 and 14 Gyr 
depending on the precise details of the HB morphology and 
the metallicity of the population. 

In general, the measured EWs of the various metal lines 
are only marginally affected by the HB extension, however 
inferred metallicities can still change if H/3 increases sig- 
nificantly, since none of the diagnostic grids is completely 
orthogonal. In practice however, these inferred changes in 
metallicity are likely to be within the observational errors in 
most cases. 



3 IDENTIFYING EXTENDED BLUE HB 
TRACERS 

One of the main issues raised by the results presented here 
is whether there are any diagnostic indices which could po- 
tentially distinguish between an old population with an ex- 
tended blue HB and a simple, single-aged population with 
an intermediate age (between, say, 3 and 8 Gyrs) in the 
integrated spectrum of an unresolved stellar system. This 
point has been raised before, but generally in the context 
of Galactic globular clusters, which all have sub-solar metal 
abundances and where it is known from isochrone fitting 
that the true ages are >10 Gyr. Here we demonstrate quite 
clearly that any old population, even at solar and super-solar 
metallicities, can have a hot HB if there is sufficient mass 
loss (coupled with a spread in mass loss) on the RGB, and 
that this will have the effect of making that population look 
spuriously young when deriving ages from Balmer lines. This 
is not merely an academic point since it is known that the 
compact elliptical galaxy M32, which has near solar metal- 
licity, has just such an extended blue H B component, which 
appe ars to centre around 9000-10000K (jBrown fc Ferguson! 
2003), the temperature at which the H/3 line is strongest. It 
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Figure 4. Fractional difference spectrum between HB model 3 
(14 Gyr plus extended blue HB) and a 3.5 Gyr SSP with r?=0.2, 
from the ct-enhanced, [Fe/H]=— 0.7 set of models. Spectra were 
normalised at 6000A. 



is not clear what fraction of the total population contributes 
to this blue HB as this is difficult to quantify from obser- 
vational data, but it is known that the spec troscopic age of 
M32, from its Balmer lines, is around 3 Gyr l|Schiavon et al.l 
120041 ). 

In order to help in the identification of potentially dis- 
criminating features in the spectra, we compared the inte- 
grated spectrum of the HB model 3 case (14 Gyr, with an 
extended blue HB, maximising H/3) with a 3.5 Gyr SSP at 
the same metallicity. We used the a-enhanced, [Fe/H] =— 0.7 
set of models for this purpose. It can be seen from Figure [2] 
that the H/3 and Fe5406 index strengths are almost exactly 
the same for the 2 scenarios, as indeed are the strengths 
of other metal indicators such as Mgfc and the various Fe 
lines (not illustrated here). The flux for the 2 spectra was 
normalised to unity at 6000A and a fractional difference cal- 
culated in the sense (3.5GyrSSP - HBmodel3)/3.5GyrSSP. 
The resulting fractional difference spectrum is shown in Fig- 
ure [4] which shows that all the differences between the 2 
cases lie at wavelengths shorter than ~4000A and, in fact, 
longwards of 4500A this difference spectrum is almost fea- 
tureless. The 3 features which obviously stand out are the 
Can feature made up of the H and K lines (3968. 5A and 
3933.7A), Mgl (2853A) and the Mgll doublet (2796A and 
2803A). It is also clear that the continuum levels shortward 
of the magnesium features are different in the sense that the 
14 Gyr population with extended blue HB has stronger UV 
continuum than the 3.5 Gyr population. 

We note here that we also performed a similar test for 
the solar metallicity, scaled solar models, comparing the 
spectrum of HB model 4 with a 3 Gyr SSP at the same 
metallicity, and found quantitatively very similar results to 
the a-enhanced case described above. 
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Figure 5. H/3 vs. Mg2800 ct-enhanced grid. The underlying grid, 
model datapoints and symbols are as for Figure [2] - only the 
[Fe/H]=— 0.7 set of models is shown, for clarity. 



3.1 The Mgn index 

At first glance it appears that the Mgn feature around 2800A 
is the most obvious candidate to investigate as an extended 
blue HB discriminator. This index , which is me a sured as an 
equivalent width, was defined by IFanelli et all (|l990l ) who 
were looking for temperature, luminosity, and metallicity 
discriminants for stellar spectra in the mid-UV. They found 
Mgn displays "no sensitivity to abundance for cool stars and 
a reversed sensitivity in FG dwarfs such that metal-poor 
stars have stronger Mgn strengths at the same temperature 
than more metal-rich stars" . We calculated the Mgn equiva- 
lent widths from the integrated spectra of all the models cre- 
ated here, as well as for the BaSTI r/=0.2 and jy=0.4 models 
used for our underlying grids. The resulting H/3 vs. Mgn a- 
enhanced diagram is shown in Figure [5] (for clarity, only the 
HB model set at [Fe/H]=-0.7 is pl otted). This g r id cle arly 
shows the behaviour referred to by IFanelli et al.l (|l990l ) as 
the higher metallicity end of the grid (i.e. towards the right 
of the diagram) folds back on itself, so that there is a large 
degree of degeneracy at old ages for all metallicities above 
[Fe/H]~ —1.0, at least for the fixed r)=0.2 case. 

More promisingly, there is a strong trend towards lower 
Mgn values as the HBs become more extended towards the 
blue. Figure [5] shows that the HB models 3 and 4 lie out- 
side the standard grid, implying a much lower metallicity 
than the true one. However, because the grid is not very or- 
thogonal, the implied age is actually close to the true age of 
14 Gyr. Interestingly the Lick index Mg6 behaves similarly 
to the Fe and [MgFe] indices displayed in Figures [2] and O 
i.e. it is negligibly affected by extended blue HBs, and so 
this index does give a good indication of the true metallic- 
ity. Hence any discrepancy between the implied metallicities 
from the Mg6 and Mgn indices for a given stellar popula- 
tion could be due to the presence of an extended blue HB 
component. 



Figure 6. H/3 vs. Call a-enhanced grid. The underlying grid, 
model datapoints and symbols are as for Figure [2] — only the 
[Fe/H]=— 1.3 set of models is shown, for clarity. Cr osses and aster- 
isks a re observational data for Galactic GCs from lSchiavon et al.l 
(2005) - the asterisks denote clusters with HB type > 0.8 (i.e. 
predominantly blue). 

There is tentative evidence that Mgn c ould be used 
in thi s way in the UV data presented in IPonder et alj 
(1998). Their Table 4 includes Mgn data for several Galac- 
tic GCs, including the well known pair M3 (NGC 5272) and 
M13 (NGC 6205), which are known to have almost identi- 
cal metallicities (around [Fe/H]=— 1.5) and similar ages, to 
within ~1 Gyr. However their HB morphologies are quite 
different, with that of M13 having a very extended blue 
tail. The IPonder et "all (|l998h data show that their Mgn 
strengths are different by 0.26 in the sense that M13 value is 
lower, which is qualitatively and quantitatively in reasonable 
agreement with our predictions. 

An important caveat to the behaviour of Mgn discussed 
here is that the current version of our population synthesis 
models does not include post-AGB stars, which are expected 
to emit much of their flux at UV wavelengths. Whilst an or- 
der of magnitude estimate indicates that their contribution 
to the continuum flux of an integrated spectrum is likely 
to be very small, it would be unwise to assume that they 
would not affect line indices. However, this magnesium fea- 
ture is clearly an intriguing possibility and merits further 
investigation. 

3.2 The Can index 

The Can index was identified by I Rose! |T984) as being highly 
sensitive to the presence of A stars in a composite spectrum. 
This index is constructed by dividing the central intensity of 
the Can H + He line by the central intensity of the Can K 
line. iRosei (I1984T ) states that "this ratio, which is constant 
in F, G, and K stars, decreases dramatically in A and B 
stars as the Can lines weaken and He strengthens". Since 
A stars have effective temperatures in the range 8000K- 
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10000K, which encompasses the range of the extended blue 
HB stars of interest here, this index merits a more detailed 
analysis. Hence we calculated Can H + He/Call K from the 
integrated spectra of all the models created here, as well as 
for the BaSTI 77=0.2 and 77=0.4 models used for our un- 
derlying grids. The resulting H/3 vs. Can a-enhanced grid 
is shown in Figure [6] (for clarity, only the HB model set at 
[Fe/H]=-1.3 is plotted). 

The first striking feature of this plot is that the lower 
metallicity 14 Gyr SSPs with 77=0.4 (the dashed line) lie out- 
side the 77=0.2 grid. Compare this with the H/3 vs. Fe5406 
grid in Figure [2] and the H/3 vs. Mgn grid in Figure [5] 
where the 77=0.4 models lie completely within the 77=0.2 
grid. Perhaps more importantly, all the extended HB mod- 
els also lie outside the 77=0.2 grid, except for HB model 1, 
which is the red clump case which has no blue component 
(and thus behaves very similarly to the fixed 77=0.2 case). 
This means that the combination of H/3 and the Can index 
could potentially distinguish between an old population with 
extended blue HB and an intermediate age or young SSP. 
Over plotted on Figu re [6] arc Galactic globular cluster data 
from lSchiavon et al.l (|2005T ) for which we have measured the 
indices directly on the observed spectra in exactly the same 
way as for our models. It is interesting to note that these 
clusters, which all have sub-solar metallicity and extended 
HBs, all lie in the region of the diagram covered by our ex- 
tended HB models. In fact, closer inspection of individual 
clusters reveals that the clusters with the bluest HBs do 
indeed fall towards the the left hand side of this diagram, 
with lower Call values. Clusters with HB type > 0.8 (i.e. 
predominantly blue) have been plotted as asterisks and it 
can be seen that they all fall well outside the fixed 77=0. 2 
'standard' grid. Looking again at the H/3 vs. Fe5406 grid 
in Figure [5] it can be seen that these same blue HB clus- 
ters would all appear to have ages around 6-7 Gyr in that 
plane. If their strong H/3 was really being caused by younger 
ages rather than blue HBs, then these clusters would have 
a higher Call ratio for the same H/3 strength and the points 
would lie further towards the right in the H/3 vs. Call grid. 

An important caveat to note regarding the Call index 
is that it is strongly affected by velocity dispersion, and also 
requires very high signal-to-noise data. This is because it 
is measured simply as the ratio of the depths at the cen- 
tral points of the Ca H and K lines, and is not an equivalent 
width (or magnitude) as most of the other diagnostic indices 
in general use are. This means it is not an ideal choice for 
studying elliptical galaxies, unless their velocity dispersions 
are known to be low. However it may be a useful indicator for 
application to extragalactic globular clusters, if sufficiently 
high signal-to-noise spectra can be obtained to measure the 
relative depths of the Ca lines accurately. Further issues re- 
lating to the Call index and its use are discussed in SectionJS] 



4 STATISTICAL FLUCTUATIONS TESTS 
4.1 Method 

In order to assess the likely impact of statistical fluctuations 
on integrated spectra, and hence diagnostic line indices, for 
real stellar systems, we simulated 2 cases representing typ- 
ical Galactic globular clusters (GCs). These 2 cases were 
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Figure 7. Diagnostic grids with results from 30 realisations of 
each of the 3 models described in the text (lines of constant age 
and [Fe/H] are as in Figure^: 2 X 10 5 M Q GC with HB model 

1 (open squares); 2 X 10 s Mq GC with HB model 3 (crosses); 

2 X 1O 4 M0 open cluster with HB model 1 (open circles). Data 
points in the lower left of each diagram are 30 realisations of the 
GC with HB model 3 modelled using 2 X 10 6 Mq, which have been 
shifted arbitrarily for clarity - as measured they are centred on 
the same location as the 2 X 1O 5 M0 data. 



based on the 14 Gyr, [Fe/H] = -1.31 (a-enhanced) SSP set 
of models listed in Table 1. Specifically we used Model 1 and 
Model 3 as the basis for these simulations, which are, respec- 
tively, the red clump HB case which does not significantly 
affect H/3, and the intermediate extension HB case, which 
results in the H/3 peak seen in Figure [5] For these tests, the 
whole of the underlying isochrone up to the TRGB, plus 
the extended HB, were populated with a discrete number 
of stars to represent a total cluster mass of 2 x 1O 5 M . For 
individual points along the RGB the mass loss was varied to 
give a spread of masses on the HB and hence a horizontal 
extension, as described in Section [2. II 

We also modelled a representative low mass cluster us- 
ing a 3.5 Gyr, [Fe/H] =— 0.7 SSP as the underlying popula- 
tion, for which we adopted a total cluster mass of 2xl0 4 M Q . 
These parameters were chosen to simulate a stellar clus- 
ter typical of those found in t he M agellanic Clouds, such 
as those studied in iDias et all (|2010T) . The reasons for this 
choice will be discussed further in Section [5] Observation- 
ally, this type of cluster would typically have a red clump 
HB (model 1) which, although extended rather than fixed 
point, is negligibly different from the fixed 77=0.2 case (as 
demonstrated in Section f2 . 2 f) . 

For all 3 of these systems, we ran Monte Carlo simula- 
tions, creating 30 realisations of each case. To explore the 
likely maximum mass at which statistical fluctuations might 
be significant in stellar clusters, we also created 30 realisa- 
tions of the 14 Gyr extended blue HB case adopting a cluster 
mass of 2 x 10 6 M Q . 
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Table 3. Icr variations in 21 Lick indices, plus 
H<5^, H7i? and Call, from the (30) extended 
blue HB simulations at 2 X 10 5 Mq. 
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A 
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A 


Fe5015 


0.106 


A 


Mgi 


0.003 


mag 


Mg 2 


0.004 


mag 


Mgb 


0.024 


A 


Fe5270 


0.054 


A 


Fe5335 


0.057 


A 


Fe5406 


0.035 


A 


Fe5709 


0.032 


A 


Fe5782 


0.017 


A 


Na D 


0.027 


A 


TiOi 


0.001 


mag 


Ti0 2 


0.001 


mag 


Can 


0.015 





4.2 Results 

Results from the statistical fluctuation tests are displayed in 
Figure [7] which shows the H/3 vs. Fe5406 and H/3 vs. Mgb di- 
agnostic grids. For the 2x 10 5 Mq GC simulations the scatter 
in H/3 is significant, and is larger for the extended blue HB 
case. However it is difficult to rigorously quantify the impact 
on implied ages for 3 reasons. Firstly, for the same degree 
of scatter in H/3, the scatter in apparent age depends on the 
absolute H/3 strength, since H/3 strength does not vary lin- 
early with age. For example, a change of ±0.2 in H/3 implies 
a shift of ~1 Gyr for a 3 Gyr SSP, but a shift of ~4 Gyr for 
a 12 Gyr SSP. Secondly, for the extended blue HB case, the 
implied age is dominated by the large offset towards younger 
ages caused by the blue HB itself, i.e. here the underlying 
14 Gyr SSP appears to have an age ~5 Gyr, even without 
the statistical fluctuations (see Figure [2}. Thirdly, the mag- 
nitude of the scatter in apparent age caused by the scatter in 
H/3 also depends on which diagnostic grid is used and what 
the properties of the underlying population are. Figure [7] il- 
lustrates this last point. In the left hand panel (i.e. the H/3 
vs. Fe5406 grid) it can be seen that the Fe5406 index also has 
a scatter of around ±0.1, in addition to the scatter in H/3. 
For the extended blue HB case this induces an uncertainty 
on both the implied age and metallicity of the population. 
However for the red HB case (model 1), the magnitude of 
the fluctuations for these 2 indices conspires to scatter the 
points along a line of constant age, implying no uncertainty 
on the age. The right hand panel of Figure [7] demonstrates 
that not all metallicity indicators behave similarly however 
- in this case, the Mgb index shows only negligible fluctua- 
tions and all the scatter is seen in the H/3 index, implying an 
uncertainty of ±2-3 Gyr for the same red HB simulations. 



As a guide to the relative significance of fluctuations in 
various diagnostic indices, la variations from the 30 sim- 
ulations of the extended blue HB case at 2 x 10 s M© arc 
tabulated in Table 3 for all 21 Lick indices, plus H<5f, H7.F 
and Call. 

As might be expected, fluctuations in H/3 are larger for 
the 2 x 1O 4 M0 case with a scatter of around ±0.25, implying 
an uncertainty of up to 2 Gyrs on the modelled 3.5 Gyr SSP. 
It should be noted that, in this case, most of the fluctuations 
are coming from discrete sampling in the upper portion of 
the RGB, since the HB consists of a very compact red clump 
with very little extension. 

Results from the 2 x 10 6 Mq simulations are also plot- 
ted in Figure [7] but have been arbitrarily shifted to the lower 
left hand corner of the diagrams for clarity. It can be seen 
that, even with this total mass, there are still fluctuations 
in the H/3 index at the level of ~0.05 in EW, corresponding 
to an uncertainty in age of 0.5-1 Gyr for old SSPs. Various 
of the metal lines also display low level fluctuations imply- 
ing a minimum uncertainty of around 0.05 dex in [Fe/H]. 
However, for these 2 x 10 6 Mq simulations all indices only 
fluctuate at a level that is likely to be within the measure- 
ment errors for real data. 



5 SUMMARY AND DISCUSSION 

To summarise, we have created integrated spectra for 16 
SSPs, 4 each at 4 different metallicities, all with an underly- 
ing age of 14 Gyr, with a range of extended HB morpholo- 
gies. This was done by varying the mass loss prescription for 
each individual SSP in 2 ways, firstly setting a mean mass, 
< M >, for stars coming on to the ZAHB, and then adding 
a spread in mass loss, om- We find that the H/3 strength 
for each SSP depends on the exact temperature distribution 
of stars along the HB, which in turn depends on the exact 
details of the mass loss prescription coupled with the metal- 
licity of the population in question. For any of the modelled 
SSPs with any amount of blue HB extension, H/3 is increased 
relative to the fixed r;=0.2 case (i.e. the 'standard' models), 
implying younger ages than the actual SSP age of 14 Gyr. 
In the worst-case scenario modelled here, a solar metallicity 
14 Gyr population with a blue HB which has a peak in its 
distribution of stars around 9000K, has an implied age of 
around 2 Gyr from the strength of the H/3 line. 

Our preliminary investigation to identify spectral fea- 
tures which might be capable of breaking the degeneracy 
between an old SSP with extended blue HB, and a truly 
intermediate a ge or young SSP, indicates that the Call in- 
dex defined bv lRosd \ 19841 ) is a very promising candidate, 
at least for populations with low velocity dispersions, such 
as extra-galactic globular clusters (but see caveat in the fol- 
lowing paragraph). There is also tentative evidence that the 
Mgn index, in combination with Mgb, could also be a very 
useful tracer of hot blue HBs. There are also indications that 
the UV continuum shortward of the Mgn feature could also 
potentially be a useful additional tool, although more work 
is needed to model this part of the spectrum in the required 
detail. However, it is interesting to note that the extreme 
blue HB models (i.e. models 3 and 4 presented here) have 
very strong UV continuum flux, as demonstrated by Fig- 
ure [4l even at solar metallicity. This supports the idea that 
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extreme blue HBs can be a significant contributor to the 
UV upturn iden tified in old elliptical galaxies (see review by 
IO'Connelllll999h . 

An important caveat to all the tests and results pre- 
sented here is that all our models have been created using 
SSPs, i.e. single age, single metallicity systems. However an- 
other scenario which could explain strong H/3 lines in a pre- 
dominantly old population would be the presence of a small 
fraction of very young stars. Although this is not likely to be 
an issue for globular clusters (Galactic or extra-galactic) it 
is a possibility which is hard to completely rule out for ellip- 
tical galaxies. Preliminary tests with our population synthe- 
sis code indicate that adding even a very small percentage 
(<1%) of a young SSP (300 Myr or less) to a 14 Gyr SSP 
would significantly strengthen H/3 and imply an interm edi- 
ate age, around 5-6 Gyr (see also ISerra fe Tragenl2007h . In 
fact, this scenario of a 'frosting' of young stars in an other- 
wise o ld stellar population has been explored bv lSmith et alj 
(2009) for (apparentl y) quiescent g alaxie s in the Shapley su- 
percluster. However. ISmith et al.l |2009l ) use the Call index 
as a tracer of hot stars and, as we have shown here, this 
index also traces hot HB stars in an SSP. So far we have 
not identified any completely unambiguous tracers that can 
distinguish between the 3 scenarios; 1) a small fraction of 
hot young stars in an otherwise old population, 2) an old 
SSP with a hot HB, 3) an intermediate age SSP. It seems 
likely that no single tracer will be able to disentangle these 3 
cases and a combination of well understood spectral indices 
and colours may well be needed - this is the focus of our 
ongoing investigation. 

Another factor potentially linked to extended blue HBs 
is the presenc e of a ste ll ar su bpopulation with a high he- 
lium fraction. iLee et all (2005) showed that the same level 
of helium enrichment required to reproduce the bluer main 
sequence in the massive GC w Centauri would also naturally 
produce the extreme blue HB stars seen in that cluster. Very 
recent observational results have provided some evidence for 
varying helium fractions within individual GCs and there is 
at least circumstantial evidence linking helium enriched sub- 
populations to extended blue HBs (see:Bra gaglia et al.ll2010l 
and references therein). The theoretical support for the HB 
morphology-enhanced-helium connection is that stars born 
with a higher helium abundance will display a lower mass 
at the turn off for a given cluster age. A lower mass at the 
TO will favour lower mass - hence hotter and bluer - HB 
stars. 

We were able to test the potential impact of a higher he- 
lium abundance on our present work by utilising the helium- 
enhanced isochrones in the BaSTI database in combination 
with a grid of helium-enhanced spectra, which have been 
produced for a subsequent paper (Coelho, Salaris & Percival, 
2011, in preparation). We chose helium-enhancement at the 
level Y — 0.3, for the [Fe/H]=— 0.7 isochrone, which is a rea- 
sonable enhancement over the standard cos mological value, 
given current observational constraints (e.g. Braaaaiia et al.l 
|2010| ) and our new spectra incorporate the same Fe and He 
abundance ratios as the isochrones. As a preliminary test, 
we measured the strength of spectral features in individual 
spectra for both the standard helium and enhanced helium 
versions, matching the spectra in T e // and loggr. We found 
that for all the diagnostic lines considered in this paper, the 
line strengths are practically indistinguishable for the two 



cases, with any differences being less than 1 percent (i.e. 
within the measurement errors). We also created an inte- 
grated spectrum for the [Fe/H]=— 0.7, 14 Gyr SSP using the 
helium-enhanced isochrone in combination with the helium- 
enhanced spectra. Again we found that the strengths of all 
the indices under scrutiny here are negligibly different from 
the standard helium c ase. We note here t hat this result is 
consistent with that of iGirardi et al.1 (|2007f ) who find that a 
similar level of helium enhancement has a negligible effect 
on broadband colours. 

The important point to stress here is that, whilst en- 
hanced helium is a likely mechanism for producing extended 
blue HBs, it in no way impacts on the work presented in this 
study. Our diagnostic models and tests are not intended to 
predict, a priori, the existence of hot HB stars, but rather 
to find diagnostics which can distinguish their presence in 
a stellar population from other hot components, such as 
a young sub-population. Whether a blue HB morphology 
arises from a large mass loss along the RGB in stars with 
'normal' helium, or from a more moderate mass loss from 
stars with higher helium, the resulting integrated spectrum 
is largely unaffected - the only parameter that matters here 
is the temperature range of stars on the HB, irrespective of 
how they have been produced. In fact, the results of our en- 
hanced helium tests, described above, show that enhanced 
helium has absolutely no impact on the diagnostic indices 
discussed here, such as H/3, Mgn and Call, nor on the results 
of the statistical fluctuations tests discussed below. 

Finally, the results of the tests presented here in Sec- 
tion 14.11 demonstrate that several key diagnostic line in- 
dices are significantly affected by statistical fluctuations, 
even in average to large mass GCs. Equally problematic is 
the fact that not all indices fluctuate at the same level - 
some are strongly affected, including all the Balmer lines, 
whilst others are only negligibly affected, such as Call (see 
Table 3). This is potentially a significant source of uncer- 
tainty if full-SED fitting methods are used to derive ages, 
metallicities and/or star formation histories from integrated 
spectra, since different features in a spectrum can be giv- 
ing confl icting best fit pa rameters. The problem is illus- 
trated in lDias et al.l (|20 10T ) who present integrated spectra 
for 14 Magellanic Cloud stellar clusters for which they de- 
rive best fit ages a nd metallicities using 2 di fferent fitting 
ades, starlight Jcid Fernandes et al.ll2005l ) and ULySS 
20091), 



<|Koleva et"aH I2009T). and 3 different sets of SSP mod- 
els, from lBruzual fe Chariot! (I2003T) .iLeBorgne et all (|2004T ) 
fPEGASE-HRl and lVazdekis et al.1 (|2010h . Half of the clus- 
ters studied have masses in the range 1 — 2.5 x 1O 4 M0, which 
corresponds to the low mass cluster modelled here, in Sec- 
tion [4TTJ For several clusters. iDias et"al] (|2010h find that the 
results from the different fitting routines and models are 
completely discrepant. As an example, cluster HW1 yields 
best-fit ages of 3.2, 5.8, 7.9, 9.0, 9.4 and 10 Gyr from the 6 
different combinations of SSPs models and SED-fitting rou- 
tines, whilst its actual age is known to be around 6 Gyr from 
isochrone fitting to the CMD. Other clusters are even more 
unconstrained, yielding ages ranging from <1 Gyr to 10 Gyr 
for the same cluster. Some of these variations in best fit pa- 
rameters are likely to be due to systematic uncertainties in 
stellar par ameters inherent within the SSP models them- 
selves fsee iPercival fe Salaris 2009), but statistical fluctua- 
tions in the observed spectra are almost certainly contribut- 
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ing to the problem. Hence we urge caution if these types of 
stellar clusters are to be used as empirical calibrating objects 
for various aspects of SPS models. 
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